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We consider a theory for a two-dimensional interacting conduction electron system with strong spin-orbit 
coupling on the interface between a topological insulator and the magnetic (ferromagnetic or antiferromagnetic) 
layer. For the ferromagnetic case we derive the Landau-Lifshitz equation, which features a contribution propor- 
tional to a fluctuation-induced electric field obtained by computing the topological (Chern-Simons) contribution 
from the vacuum polarization. We also show that fermionic quantum fluctuations reduce the critical temperature 
T c at the interface relative to the critical temperature T c of the bulk, so that in the interval T c <T < T c is possible 
to have coexistence of gapless Dirac fermions at the interface with a ferromagnetically ordered layer. For the 
case of an antiferromagnetic layer on a topological insulator substrate, we show that a second-order quantum 
phase transition occurs at the interface, and compute the corresponding critical exponents. In particular, we 
show that the electrons at the interface acquire an anomalous dimension at criticality. The critical behavior of 
the Neel order parameter is anisotropic and features large anomalous dimensions for both the longitudinal and 
transversal fluctuations. 
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A spin current may exhibit interesting topological proper- 
ties in systems where a Berry curvature in Bloch momentum 
space is induced by the underlying band structure Q]|2), ^ e 
for example in the case of some hole-doped semi-conductors 
described by a Luttinger Hamiltonian O or systems featur- 
ing a Rashba spin-orbit coupling Q. More recent prominent 
examples involve the quantum spin Hall insulators or topolog- 
ical insulators (TI) J5]|6l, where a Berry curvature in momen- 
tum space also arises. Depending on the physical situation the 
Berry curvature may be Abelian or non-Abelian, and deter- 
mines a magnetic monopole in momentum space. 

The surface of a TI, when in contact with a material ex- 
hibiting magnetic order, offers a framework for many topolog- 
ical effects. For instance, the two-dimensional system repre- 
sented by the surface of a topological insulator can be used as 
substrate for a magnetic layer, which can be either ferromag- 
netic (FM) or antiferromagnetic (AF). For a FM layer hav- 
ing a TI as substrate, a theoretical study of the magnetization 
dynamics was carried out recently [7||. In a similar context, 
the electric charging of magnetic textures has also been dis- 
cussed [8]. Other interesting electromagnetic topological ef- 
fects with a similar setup were studied lr9 UT2l and have shown 
to exhibit properties similar to those of axion electrodynamics 
ITOI . In the axion electrodynamics a topological term of the 
form (%n 2 y x a6R ■ B is present IfTUl [131 in the action, where 
a is the fine structure constant and 9 is the so called axion 
field. For the case where 9 is uniform, time-reversal invari- 
ant TIs require 9 — n ifTOll . Such a term should play a very 
important role at interfaces of TIs with other insulators. For 
a magnetic insulating layer on the surface of a TI, a modifi- 
cation of the magnetization dynamics occurs, due to a direct 
coupling of the magnetization to the electric field. Indeed, we 
have E • B = E ■ (H + AnM), giving rise to a magnetoelectric 
effect, which influences the precession of the magnetization 

At the same time, the experimental situation is far from be- 



ing clear, as many theoretical predictions have not been con- 
firmed yet. Moreover, in some cases, it seems that the theoret- 
ical expectations are not fulfilled. For instance, from a theory 
perspective one would expect that the coupling of a TI to a 
FM layer makes the surface states gapped. However, in a very 
recent experiment [ 14 1 where Fe impurities were deposited on 
Bi2Se3, no sign for a gap has been found, in apparent conflict 
with theoretical expectations. Therefore, further theoretical 
studies on the coupling of a TI substrate to magnetic system 
are necessary. 

In this paper we consider quantum fluctuations effects stem- 
ming from the proximity-induced magnetism on the surface of 
a TI. We assume that the electrons on the surface of the TI in- 
teract via a long-range Coulomb interaction. For the case of 
a FM layer in contact with the TI, we will derive a Landau- 
Lifshitz (LL) equation which accounts for these interaction 
effects. In our calculation an axion-like term emerges due to 
quantum fluctuations. At the interface, it manifests itself as 
a Chern-Simons (CS) term 1151 . which breaks time-reversal 
symmetry, as a consequence of the coupling of the surface 
of the TI to the magnetic layer. Furthermore, the electronic 
quantum fluctuations make the stiffness anisotropic, even if 
the bulk of the FM layer features an isotropic stiffness. We 
also show that due to the quantum fluctuations of the elec- 
trons, the critical temperature T c at the interface is reduced 
relative to the critical temperature T c of the FM layer. This 
allows the existence of gapless fermions at the interface at the 
temperature range f c < T < T c where the bulk magnetic layer 
is still magnetically ordered. 

It has been shown in a recent study (16 1 that the best can- 
didate material to gap the topological surface states of a TI 
is MnSe, which is an AF insulator. Therefore, the study of 
proximity-induced AF order on a TI is also of great relevance. 
For such a case we will show that a second-order quantum 
phase transition occurs at the interface, and that it defines 
a new universality class. One consequence of this interface 
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quantum criticality is that the surface electrons become gap- 
less at the quantum critical point (QCP). This does not hap- 
pen in the FM case we study, where the electrons are always 
gapped at zero temperature. Hence, the topologically pro- 
tected gapless modes can be restored at zero temperature by 
disordering the AF long-range order at the interface. A fur- 
ther important feature of this interface quantum criticality is 
the emergence of a large anomalous dimension for the Neel 
order parameter. Interestingly, at the QCP the fermions will 
also be shown an anomalous dimension. Large anomalous di- 
mensions are known to occur in systems exhibiting phase tran- 
sitions that do not fit the Landau-Ginzburg- Wilson paradigm, 
like in the case of the so called deconfined quantum criticality 
scenario in some Mott insulators ifTTl . In this context, large 
anomalous dimensions in theories involving Dirac fermions 
have also been discussed [ 18l [T9l . 

Our starting point is the Lagrangian for conduction elec- 
trons interacting via a Coulomb interaction on the surface of 
an insulator either in contact with a bulk FM composed of sev- 
eral layers, similarly to Ref. 0, or with an AF bulk system. 
Thus, if n is the induced magnetization at the interface and 
L the angular momentum, the spin of the conduction elec- 
trons, S = (l/2)c' crc, is coupled to the total magnetization 
(jUb/2)L + n via an exchange term -2/oS • [(//b/2)L + n/Ay, 
where = [cj cj], with & - (o- x ,cr v ,cr z ) being the Pauli 
matrices, and Ni is the number of magnetic insulating layers. 
The lack of inversion symmetry in the direction perpendicular 
to the interface leads to a spin-orbit coupling of the Rashba 
type. Thus, the Lagrangian for the conduction electrons at the 
interface reads (we are assuming units where ti — c — 1), 



£, c — c [idt + fi + eip — iv(p~ y d x - cr x d y ) + Jcr ■ n]c 



(1) 



where v oc J and J = Jq/N^. In writing the above Lagrangian 
we have assumed that the spin-orbit coupling is much stronger 
than the usual kinetic term of the conduction electrons, which 
has been neglected. The auxiliary (Hubbard-Stratonovich) 
field <p accounts for the Coulomb interaction. Upon integrat- 
ing out <p the usual Coulomb interaction between the electrons 
is obtained. The non-local Gaussian term for <p reflects the 
three-dimensional character of the Coulomb interaction in a 
two-dimensional problem, similarly to graphene ll20l . In this 
term e represents the dielectric constant. 

The full Lagrangian of the systems includes the Lagrangian 
describing the magnetization dynamics, which for a FM sys- 
tem includes a Landau-Ginzburg (LG) type functional and is 
given by 



£ FM = b-d,n-^(Vn) 2 -^n 2 - J(n 2 ) 2 , 



(2) 



where k, u > and m 2 = ao(T - To), with To being the 
(mean-field) critical temperature to disorder the FM. b is the 
Berry connection, which fulfills the usual monopole condi- 
tion, dbj/drtj - dbj/drii = e^nul^- For m 2 < (or T < To) 
the bulk FM is in a ferromagnetically ordered state. 



Before considering the magnetization dynamics, let us first 
consider a fluctuation-corrected mean-field theory where cr is 
assumed to be uniform, and that the transverse fluctuations 
vanish. We also neglect the fluctuation effects of the Coulomb 
interaction. The calculations are done in imaginary time and 
finite temperature. 

In this case, after integrating out the fermions, we obtain 
the free energy density, 



T 



^ I C 2^ ln(w " + v2p2 + ^ ) + Y cr " 



-tY 



+ — cr, 

4! 

(3) 

where u>„ = (2m + \)jiT is the fermionic Matsubara frequency. 
After performing the Matsubara sum, the remaining integral 
over momenta contains a zero temperature contribution which 
is divergent, requiring regularization and renormalization. Us- 
ing an ultraviolet cutoff A ~ aT x , where a is the lattice con- 
stant, we can cancel the dependence on the cutoff by mini- 
mally absorbing it in a redefinition of the Curie temperature of 
the bulk FM precisely at the interface. The physical require- 
ment (or renormalization condition) is that the zero tempera- 
ture fermionic gap, = 7cr , is finite in the long-wavelength 
limit. 

The saddle-point approximation yields, 



u , J l T 
ao(T c -T)= -cr z + — - In 



2 cosh 



(27) 



(4) 



where T c is the renormalized Curie temperature of the bulk 
FM at the interface. The critical temperature, f c , at the in- 
terface is obtained by demanding that cr vanishes at T — T c . 
This yields f c = T c [l + J 2 lnl/C^oV 2 )] -1 . On the other hand, 
by setting T = in Eq. Q, we obtain that at the interface 
T c = um 2 /(6J 2 ) + J 2 m,/,/(27raov 2 ). Note that this expression 
makes only sense at the interface and does not correspond 
to the physical critical temperature there, which is actually 
given by T c , Furthermore, since v oc J and at leading order 
cr 2 » 6aoTo/u, we obtain that T c — > To as v — > 0. In order 
to estimate T c , we assume that T c » um 2 /(6J 2 ao), such that 
we have approximately, T c * m^,T c [m^ + (21n2)7 , c ]~ 1 . If we 
use the estimates m 28.2 meV and T c » 70 K irT6ll2TI . we 
obtain T c =s 54 K. Thus, our fluctuation-corrected mean-field 
theory implies that the critical temperature at the interface is 
smaller than the Curie temperature of the bulk FM. Therefore, 
it is possible to destroy the proximity-induced magnetization 
at the interface while the bulk FM is still ordered. This oc- 
curs typically in a temperature window T c < T < T c , where 
we are assuming that T c does not differ appreciably from To- 
The reduction of the critical temperature at the interface with 
respect to the bulk one is an important consequence of the in- 
terplay between ferromagnetic proximity effect and the spin- 
orbit coupling. 

The next step is to compute the fluctuations of the order 
parameter around the mean-field theory. Since we are inter- 
ested in deriving a differential equation for the magnetiza- 
tion dynamics, we will return to real time in the following 
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FIG. 1. Diagrams contributing to 6S gmge , Eq. (jTJ. The wiggled line 
represents the vector field a^, the solid line is a Dirac fermion, and 
the dashed line represents the fluctuating part of the cr field. The 
diagrams (b) and (c) cancel out. 



and essentially consider a zero temperature calculation. In or- 
der to facilitate our analysis of the problem, it is convenient to 
rewrite the Lagrangian for the conduction electrons in a QED- 
like form, which is achieved by the rescalings if — > (J/e)<p, 
Xj — > vxj (i = 1,2), in the action, to obtain 

Lc = tA0$ - JM + Jcn/fift - -flo|V|flo, (5) 

where the Dirac matrices are defined by y° = cr z , y l = — icr x , 
and y 2 = icr y , i[> = vc, such that the usual relativistic notations 
for spinors hold with \jj = ifry , and also the usual Dirac slash 
notation, (g = y^Qf,, is being used. The gauge field is given 
by cf = (tp, n y , -n x ) and cr = n z , and the dielectric constant, 
f = ev 2 J 2 l(2jie 2 ). The chemical potential was absorbed in a . 
Note the different physical nature of the components of <f, 
with the zeroth component corresponding to the scalar poten- 
tial, and the spatial components being related to the transverse 
components of the magnetization relative to the z-direction. 

The main advantage of the QED-like representation of the 
Lagrangian for the conduction electrons is that it allows us 
to integrate out the fermions using standard field theoretic 
techniques. In order to better explore different approximation 
schemes, we will generalize our original Lagrangian to the 
case where there are fermionic orbital degrees of freedom. 
Thus, integrating out the fermions yields the gauge-invariant 
contribution to the effective action, 



iNTr\n(ifi- J<ji + Jo-). 



(6) 



The lowest order diagrams associated to the fluctuating 
fields from Eq. |6]l are shown in Fig. [T] Approximate evalu- 
ation of S g aU ge at long wavelengths yields the leading fluctua- 



tion contribution, S r 



6S, 



NJ 2 



;e ~ + 55 gauge , where 



+ e^SV + — [(d,&Y - (Vcr) 1 ] 



(7) 



with the fluctuation <x arising from the decomposition cr = 
(cr) + <x. The quadratic fluctuation term in & will generate an 
anisotropy in the magnetic system, which is isotropic in the 
bulk. The first term in Eq. |7]) corresponds to a Maxwell term 
in (2 + l)-dimensional electrodynamics. The second term is 
a CS term [15] generated by the quantum fluctuations. This 



CS term reflects the breaking of time-reversal symmetry due 
to the coupling to a magnetic layer. In order to better appre- 
ciate the effect of the CS term, it is useful to rewrite the CS 



contribution to S gauge in the form, 



> cs 



An 



! 



d i x(n y d t n x — n x d,n y + 2n ■ Vcj), (8) 



where cr xy = cr° xy NJ 2 /e 2 , with cr° Ky = e 2 /2 (in units where 
ft = 1), is the induced Hall conductivity. It is readily seen 
that the contribution proportional to n x d,n y - n y d t n x yields an 
additional Berry phase, as discussed previously in Ref. J7). 
The term proportional to n ■ Vip is a crucial contribution stem- 
ming from the Coulomb interaction between the fermions at 
the interface. Indeed, since if is a fluctuating scalar poten- 
tial associated to the Coulomb interaction, this term yields 
a contribution proportional to M • E, where E = -Vif is a 
fluctuation-induced electric field, and M ~ n. Thus, this term 
corresponds to an emergent axion-like term. 

From the CS term we can also determine the topological 
current, = (4n)~ l cr xy e MvA d v aA, which up to the prefactor 
<j XY l{An) has components (V ■ n,d t n x - d y <f,d,n x + d x <f). The 
topological current is identically conserved, i.e., d^W = 0. 

The Landau-Lifshitz equation for the magnetization dy- 
namics at the interface can now be obtained by extremizing 
the effective action (note that here V • n = d x n x + d y n y , since 
there is no z-component for the V operator in the present two- 
dimensional problem), 

d,n = n x (p'/CV 2 !!)^,- + ).c- - ViV ■ n)| 



Zcr xy 

2KV 2 



n xE 



3«v 



(n ■ e,)<9 f E 



(9) 



where the stiffness matrix elements are given by p'J = 
(Z/v 2 )[K(6 ix 6 jx + 5 iy 5jy) + {k + cr^/ilnm^di.dj,], with Z = 
[1 - m^(T xy l{2nv 2 J)y x . Eq. (|9jl is one of the main results 
of this paper. It leads to a fluctuation-induced magnetoelec- 
tric effect. One important consequence of Eq. (|9| is that due 
to the fluctuation-induced electric field, the magnitude of the 
magnetization is not constant, as it would be in the case of ab- 
sence of Coulomb interaction or for a constant electric field. 
In particular, if the electric field is only due to external effects, 
this result implies that we can use a time-dependent electric 
field to control the magnitude of the magnetization. 

Part of the coupling to the electric field, discussed previ- 
ously by Garate and Franz 0, is reproduced here as a fluctu- 
ation effect due to the Coulomb interaction between the spin- 
orbit coupled electrons lying on the surface of a TI. We have 
obtained in addition a contribution involving d,E that accounts 
for the time-dependence of the electric field. Note that the 
term involving (d 2 n) z is typically small at low energy and can 
be safely neglected in most calculations. 

Next we discuss the case of an AF layer on a TI substrate 
at zero temperature, which, as we will see, differs fundamen- 
tally from the case of a FM layer. In the AF case a quantum 
phase transition occurs at the interface. In order to study the 
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phase structure of the theory in this case, we will work only in 
imaginary time from now on. Specifically, we consider an Eu- 
clidean effective field theory whose Lagrangian has the form, 

/ 1 , , 

£. = i//($- igi4 + 82(r)iff + -fl |V|fl + ^ [(^cr) + (df,a) ] 



M 2 2 2\ ^ / 2 1 1 

+ ~y(o- z + a 2 ) + —{a 1 + a 2 ) 2 , 



(10) 



where we are not making any longer distinction between up- 
per and lower covariant indices, since the metric of the theory 
has now Euclidean signature. Note that we are also assuming 
that g\ + g2, as quantum fluctuations induces an anisotropy. 
In the spirit of effective field theories, the coupling constants 
are understood as effective parameters to be determined by a 
renormalization group (RG) flow. Thus, the phase structure of 
the theory is completely determined by the RG equations for 
the coupling constants. 

At low energies and one-loop order (see below) the fixed 
point structure will be governed by the dimensionless cou- 
plings g 2 = g 2 r /M 6 r (i = 1,2) and A = A r /M e r , where A r and 
gi r are corresponding renormalized couplings and we are us- 
ing the renormalized mass M r as renormalization scale El . 
Here e = 4 - d, where d = D + 1, with D being the spatial 
dimension. Our analysis is done in the framework of the e- 
expansion, which is carried out up to one-loop order. As usual, 
in such a renormalization scheme the renormalized mass gives 
the inverse of the correlation length, i.e., M r = Due to the 
coupling between <x and the fermions, a mass anisotropy will 
be generated, defining in this way two correlation lengths, re- 
lated to longitudinal and transversal fluctuations. We will as- 
sume that ^ refers to the longitudinal correlation length, giv- 
ing the fluctuations of the <x field. The correlation length 
due to transversal fluctuations will be denoted by g x . If v 
and v± are respectively the critical exponents of the longitu- 
dinal and transversal correlation lengths, we easily obtain that 
£ ~ £ v±/,v , which determines the crossover exponent <p = v/v±. 
The quantum critical behavior can be derived from a gener- 
alization of the extended Gross-Neveu model |23| discussed 
in Ref. [24 1. We obtain in this way the one-loop p functions 
P~ g i = M r dg 2 JdM r , Pg = M r dg\ldM r and f3 x = M r dA/dM r in 
the form, 



.2 N 
+ 12^ 



N + 3 



8tt 2 



Pa 



-eu + 



1 /H 



8^ 2 



■2NAgl 



UNg 4 2 



(11) 



(12) 



(13) 



The p function for ( = ( r /M r follows from the non-locality 
of the quadratic term in «o- Since counterterms are local, 
this term does not renormalize, which implies simply Pf = 
(Ng 2 /I2n 2 - e)l. Note that at one-loop order the RG flow for 
g ? decouples from the other couplings. 



The quantum critical point is determined by demanding that 
the p functions above vanish, which yields the infrared stable 
fixed points, g 2 u = I2n 2 e/N, g 2 t = 8n 2 e/(N + 3), and X = 
8^ 2 e(3 -N + VA^ 2 + 258JV + 9) /[ll(iV+3)]. The anomalous 
dimension rj^ of the Neel order parameter at the interface can 
be defined via the scaling behavior (<x) ~ M (2 e+nN)/2 ^ anc [ j s 
given at one-loop by rj N = Ng 2 lt l(%n 2 ) = Ne/(N + 3). For 
N = 1 and two spatial dimensions (corresponding to e = 1), 
we obtain rj N = 1/4. This large value of the anomalous di- 
mension, as compared to the value obtained from the 0(3) 
universality class, reflects the fact that {&) receives contribu- 
tions from the composite operator i/n/'. The scaling behav- 
ior of n = (n x , riy, a) at the interface is anisotropic, and the 
transversal fluctuations have a different anomalous dimension, 
which is dominantly determined by the vacuum polarization 
diagrams, rj^ = e, yielding rjj, = 1 for D = 2. It is worth to 
mention that two-loop corrections will be small, but positive 
(typically ~ 0.03). Therefore, we expect that a more accurate 
value for the anomalous dimension ?^ is slightly above the 
unity. 

The electrons at the interface also have an anomalous scal- 
ing at the quantum critical point, where (cr) vanishes at zero 
temperature and the electrons become gapless. This is in 
contrast with the FM case, where the fermionic spectrum 
is always gapped at zero temperature. Thus, we obtain the 
low-energy behavior, (tff(p)t//(p)) ~ -ip/p 2 ^*, where = 
g 2 2 J(l67r 2 ) = e/[2(N + 3)]. For D - 2 and N = 1, we obtain 
r\ty — 1/8. Note that 77^ does not receive any contribution from 
the fixed point g 2 1 at one loop order. This is due to the fact 
that the vector field propagator takes here the same form as 
one in QED where the Feynman gauge has been fixed. 

It remains to compute the critical exponents of the cor- 
relation lengths. The longitudinal correlation length expo- 
nent is given by v = (2 + t]m) 1 , where at one-loop, t]m = 
-5/l*/(487r 2 ) - t]n. Thus, by expanding up to first order in e, 
we obtain, 



1 e 

2 + 4(JV + 



3) 



5 

66 



(3-N+ y/N 2 + 258JV + 9j + N 



(14) 

Setting once more D — 2 and N - 1, we obtain v as 0.649. The 
transversal correlation length exponent, on the other hand, is 
given by v ± = (2 + 77^), where rfa = t] M ~ In + Iff- Thus, 
we obtain v± ~ v + 3e/[4(N + 3)], which for N - 1 and 
D — 2 yields v± ~ 0.83. Note that the values of the corre- 
lation length exponents differ appreciably from the one-loop 



value of the 0(3) universality class, 



one-loop 

v OO) 



0.614. The 



correlation length exponents calculated at one-loop are con- 
siderably less accurate than the anomalous dimensions, since 
two-loop corrections usually make significant changes in v. 
However, our calculation clearly shows the trend towards a 
novel quantum universality class at the interface of a TI with 
an AF. 

In conclusion, we have shown that the interface between a 
topological insulator and a magnetic layer exhibits interesting 
fluctuation-induced magnetoelectric effects. Indeed, we have 
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shown in the FM case that an axion-like term is generated in 
the form of a CS term, which in turn modifies the magneti- 
zation dynamics of the LL equation. Furthermore, we have 
shown that for a specific temperature window is possible to 
have gapless fermions at the interface and, at the same time, 
a ferromagnetically ordered layer. It would be interesting to 
check this result experimentally by studying the temperature 
variation of the surface states. 

For the case of an AF layer, we have shown that a quantum 
phase transition occurs at the interface, and that the fermion 
spectrum becomes gapless at the QCP. Moreover, large val- 
ues of the anomalous dimensions for the Neel order parameter 
were obtained. 
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